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ABSTRACT 

A long 280 ks observation of the Seyfert 1 galaxy NGC 3783 with XMM-Newton 
is reported. We focus on the oxygen line complex between 17 and 24 A as mea- 
sured with the RGS spectrometers. Accurate absorption column densities and 
emission line fluxes are obtained. We explore several options for the geometry 
and physical form of the emitting and absorbing gas. The lack of change in ioniza- 
tion in the absorber despite an increase in continuum flux during the observation 
restricts the high-ionization (O-K) and the low-ionization (Fe-M) gas to distances 
of at least 0.5 pc and 2.8 pc, respectively, away from the central source. Given 
the P-Cygni type profiles in the resonance spectral lines and the similar velocity 
widths, column densities, and ionization structure inferred separately from the 
emission and absorption lines, it is tempting to relate the X-ray narrow-line emit- 
ting plasma with the X-ray absorbing gas. Under this assumption, the scenario 
of dense clumped clouds can be ruled out. Conversely, extended ionization cones 
(r > 10 pc) are consistent with the observation independent of this assumption. 
These findings are in stark contrast with the picture of numerous clumpy (n e > 
10 9 cm -3 ) clouds drawn recently from UV spectra, but it is consistent with the 
extended X-ray emission cones observed directly in Seyfert 2 galaxies. 
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1. Introduction 

As a bright, nearby Seyfert 1 galaxy [z = 0.00976, de Vaucouleurs et al. (1991)], 
NGC 3783 is one of the best studied Active Galactic Nuclei (AGN). In the X-rays, observations 
with ROSAT (Turner et al. 1993) and ASCA (George et al. 1998) revealed a highly-ionized 
absorber, but lacked the spectral resolution required for kinematic measurements. Lately, 
the rich absorption line spectrum of NGC 3783 has been readily measured with the grating 
spectrometers on board the X-ray observatories Chandra (Kaspi et al. 2000, 2001, 2002) and 
XMM-Newton (Blustin et al. 2002), telling us that the highly- ionized gas is outflowing at 
velocities of a few 100 km s -1 . These spectra of NGC 3783 feature numerous absorption 
lines by K-shell ions of C, N, O, Ne, Mg, Si, and S, as well as by L-shell ions of Fe. Many 
inner-shell absorption lines formed by lower ionization states are also found in NGC 3783 
(Behar & Netzer 2002) indicating a broad range of ionization. 

Kaspi et al. (2002) used the high resolving power (A/AA up to 1000) of the Chandra 
High Energy Transmission Grating Spectrometer (HETGS which covers the 1.6 - 23.4 A 
wavelength range) to determine the mean outflow velocity of —590 ± 150 km s -1 and a 
mean FWHM of 820 ± 280 km s _1 for the absorber in NGC 3783. The unprecedented small 
errors are due to an integrated observing time of 900 ks. Both velocities are consistent with 
a blend of known UV absorption systems. Blustin et al. (2002) employed a short 40 ks 
observation of NGC 3783 with the Reflection Grating Spectrometers (RGS, A/AA up to 400) 
on board XMM-Newton to construct a plasma model with two ionization phases in the X-ray 
absorber that reproduced the RGS spectrum between 6 and 38 A quite well. However, as 
was also pointed out by Kaspi et al. (2001), the wide range of observed charge states is 
indicative of a broad distribution of ionization within each phase. Despite the extensive X- 
ray measurements, the location, form, and energetic significance of the outflow in NGC 3783 
have remained much of a mystery. 

The connection of the X-ray absorber to AGN constituents observed in other wavebands 
is a long standing puzzle. The broad-line region (BLR) in NGC 3783 is inferred from the 
optical time lag measurements of Onken & Peterson (2002) to cover a range of distances 
from approximately 10 15 to 10 16 cm (not along the line of sight). Combining X-ray and UV 
observations, Shields & Hamann (1997) favor the interpretation by which the UV and X-ray 
absorbers constitute a single medium. Their analysis implies an upper limit of about 30 pc 
for the separation between the nucleus and the absorbing plasma. Recently, four velocity 
components have been resolved by Gabel et al. (2003) in the UV absorption spectrum of 
NGC 3783. Interpreting absorption of C III at 1175 A as arising from a metastable level, 
Gabel et al. (2003) estimate the electron density in the high-velocity component to be n e = 
10 9 cm~ 3 and, hence, the absorber to be located at a distance r < 8 x 10 17 cm away from 
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the central source and with a typical size of Ar = 10 10 cm, i.e., numerous highly-clumped 
cloudlets. Conversely, a narrow O VI emission line component with a centroid velocity of 
km s -1 was estimated by Gabel et al. (2003) to originate about 40 pc away from the 
nucleus. 

In December 2001, XMM-Newton observed NGC 3783 for 280 ks. The high effective area 
of the RGS makes this data set comparable in statistical quality to that obtained in 900 ks 
with the HETGS in the respective wavebands. In this paper, we present the RGS spectrum 
from this long observation. The main goal is to further characterize the X-ray absorbing and 
emitting gas and particularly to constrain its geometry and physical conditions. 

2. The Data 

2.1. Data Reduction 

NGC 3783 was observed by XMM-Newton for a total of 280 ks from December 17 - 21 
as part of the RGS guaranteed time program. The data obtained with the EPIC X-ray 
cameras and with the Optical Monitor (OM) will be discussed elsewhere (Blustin et al., in 
preparation). Here we focus on the RGS spectrum. The RGS data were processed using 
custom software developed at Columbia University, which operates on the observation data 
files (ODF). This software is similar in origin and in function to the RGS branch of the 
science analysis system (SAS), but is more flexible and includes recent calibration of the 
instrumental O-edge region based on accumulated data from long observations of Mrk 421 
(Rasmussen, Paerels, & Kahn 2003). It corrects for the effect of aspect drift on the dispersion 
angles and includes background subtraction using off-target CCD regions. 

2.2. Short Time Variability 

The RGS light curves are shown in Fig. 1. It is evident that the soft X-ray flux varies 
fitfully and by as much as a factor of 2 on time scales of hours. The light curves obtained 
with the EPIC cameras, which are more sensitive to hard X-rays, show very similar behavior. 
The UV light curve obtained with the Optical Monitor (OM) resembles the curve in Fig. 1, 
only it is somewhat smoother (Blustin et al., in preparation). Observed variations in the 
ionization state of the absorber due to changes in the continuum flux could potentially shed 
light on the location of the absorber through ionization and recombination times. With 
that goal in mind, we divide the 280 ks observation into four temporal segments, seeking 
variations in the absorption spectrum. 



-4- 



The four segments are labeled in Fig. 1 as the "low", "rise", "high", and "peak" phases. 
The spectra obtained for each segment separately are presented in Fig. 2. Each segment 
by itself suffers from limited statistics. Consequently, disparities in narrow features such 
as emission or absorption lines, which have limited counts in a few bins, are very difficult 
to detect unambiguously. On the other hand, it is clearly seen that the broad absorption 
features in the spectrum, such as the Fe M-shell 2p - 3d unresolved transition array (UTA) 
between 16 and 17 A and the O VII and O VIII photoelectric bound-free K-edges (just 
noticeable at 14.2 & 16.8 A) do not vary significantly from one time segment to the other. 
Indeed, calculated ratios of spectra among the various time segments show mostly noise 
around the flat ratio of the continua to within the data errors, i.e., ~10%. For example, the 
ratio of the "high" and "low" segments is plotted in Fig. 3. The plot is rather noisy, but 
no clear bound-free edge features are observed. The only statistically significant deviations 
from the mean in Fig. 3 are well localized around strong emission/absorption lines (e.g., 
Ne IX at 13.48 A). Since emission lines are not expected to disappear in the ratio plot and 
since ratios at deep absorption troughs are unreliable, we interpret the lack of bound-free 
features as evidence for no detectable variation in absorption. 

The position and shape of the Fe-M 2p-3d UTA are expected to be particularly sensitive 
to the ionizing flux as several M-shell charge states of the same element are observed at the 
same time. During the present observation, the ionizing flux overall increased gradually by 
almost a factor of 2 (Fig. 1). An average increase of about 50% occurs from the "low" state 
to the "high" state. Such an increase in the ionization parameter £ [£ = L/(n e r 2 ), where L is 
the ionizing luminosity] would have shifted the UTA towards shorter wavelengths by about 
0.4 A, if the gas had reached a new equilibrium state (c.f., Fig. 4 in Behar, Sako, & Kahn 
2001). Such a shift would have been surely detectable with the RGS. Carefully examining 
Fig. 3, it is possible to visualize a slight, statistically insignificant, shift in the UTA around 
16 A. However, this is exactly where the O VIII Ly/3 line falls and in any event the effect 
is smaller than would be expected from new equilibrium conditions. Together with the lack 
of edge features in Fig. 3, this makes for strong evidence that the ionization state in the 
absorber does not have ample time to reach a new equilibrium state and, thus, it is hardly 
affected by the rapidly varying continuum. 

Since the ionization time depends directly on the absolute ionizing flux seen by the 
absorber (i.e., its distance from the source), we can use the lack of ionization despite the 
increase in flux to place a lower bound to the separation between the source and the ab- 
sorber. For a detailed discussion of ionization in winds, see Krolik & Kriss (1995). We 
take the average power-law ionizing flux density F E = A(E/E )- r with E = 1 keV and 
r = 1.7 (determined for NGC 3783 by Kaspi et al. 2001; Blustin et al. 2002, and here). 
The normalization factor A scales with r~ 2 and reaches its observed value of about 0.015 
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photons s^ 1 cm" 2 keV -1 at the telescope at a distance of 41.8 Mpc (assuming H = 70 
km s~ x Mpc^ 1 ). To first order and as an upper limit, the equilibration time is the reciprocal 
of the photoionization (PI) rate J Fe<j pi (E)dE, where a PI (E) is the PI cross section. For 
O VII, we get a time of 4.6 r? c , days, where is the distance of the absorber from the 
source in parsecs. The unnoticeable change in ionization observed over more than a day 
following the rise of the continuum implies that the O VII absorber is at least 0.46 pc away 
from the source. Extrapolating the power law continuum down to the ionization energy of 
Fe IX (= 0.23 keV), which is dominant in the Fe-M UTA (at about 16.6 A), but also repre- 
sentative of the other M-shell ions, we get an even shorter ionization time of about 0.13 r 2 ^ 
days. Consequently, the Fe IX absorber must be at least 2.8 pc away from the source to 
not be affected. Although only crude lower limits, these estimates already clearly exclude 
the BLR, which is approximately 0.003 pc away from the source (Onken & Peterson 2002). 
For comparison with previous estimates, Appenzeller & Wagner (1991) place emission by 
highly-ionized species (e.g., Fe VII) just outside the BLR, while Shields & Hamann (1997) 
give an upper limit of 30 pc for the X-ray and UV gas. 

Note that despite the fluctuations over time scales of days, the average X-ray flux in 
NGC 3783 has remained steady over the past few years (e.g., Markowitz et al. 2003). Under 
these conditions and at a distance > 0.5 pc, the ionization state of the absorber does not 
react to the short time-scale fluctuations, but instead it is governed entirely by the mean 
flux. Indeed, Kaspi et al. (2001) point out that the shape of the X-ray absorption spectrum 
of NGC 3783 is unchanged since 1996 when it was measured with ASCA, although the 
continuum level has clearly varied during that period. A study of spectral variability on 
intermediate time scales of a few months using fractions of the 900 ks HETGS observation is 
bound to shed new light on the behavior of NGC 3783 in time (Netzer et al., in preparation). 



2.3. Time Integrated Spectrum 

The fact that the absorption in the spectrum does not vary significantly during the 
course of our observation justifies the treatment of the 280 ks time integrated spectrum as 
a whole for studying the absorption. Furthermore, when we compare the present average 
spectrum with earlier observations (Kaspi et al. 2001, 2002; Blustin et al. 2002), we still find 
no significant discrepancies. Consequently, for the most part the best-fit model proposed 
by Blustin et al. (2002) describes the current spectrum rather well. This is demonstrated 
in Fig. 4 where the slab model from Blustin et al. (2002) is compared with the current 
280 ks spectrum. No adjustment whatsoever has been applied to the model, including no re- 
normalization. Noticeable discrepancies between data and model occur in the details around 
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the O Ka region between 14 and 24 A. These could not have been modeled properly with 
the first 40 ks exposure due to lack of counts in the narrow lines present in this range. 

In this work, we focus on the O I - O VIII complex in the wavelength range 17.3 - 
24.7 A. The goal is to measure individual emission and absorption line fluxes. Due to 
considerable blending, a simultaneous fit of many spectral features is required. We fit the 
spectrum locally with a power-law plus an ensemble of individual-ion absorption models and 
a few (unabsorbed) emission lines. Subsequently, galactic absorption (Nh = 8.7xl0 20 cm" 2 ) 
is applied to the entire spectrum. The models are based on atomic data calculated with 
HULLAC (Bar-Shalom, Klapisch, & Oreg 2001) and include all of the inner-shell atomic 
transitions and their autoionization widths for all ions and elements that appear in the 
spectrum. A similar method has been successfully applied to several previous AGN spectra 
(e.g., Sako et al. 2001; Kaastra et al. 2002; Blustin et al. 2002). 



2.3.1. Absorption 

Out of four velocity components resolved in the UV (e.g., for O VI), of which two are 
marginally resolved with HETGS in lines of O VII and O VIII (Kaspi et al. 2002), the RGS 
reveals only one blended trough. Indeed, the blueshifts measured here (= —600 km s _1 ) are 
consistent with the mean velocity of the two dominant (& slowest) components in the UV 
and are in agreement with the velocity measured with HETGS. Thus, in the model, we use 
only one velocity component for each ion. We fit for the turbulent velocity width and obtain 
Vturb — 170 km s _1 (i.e., o~ v = 120 km s _1 , FWHM = 280 km s _1 ). This width is very close 
to the value of 150 km s _1 used by Kaspi et al. (2000). Both here and in Kaspi et al. (2000), 
this velocity width was used because it enables the fitting of entire series of lines with a 
single kinematic component. However, it should be clear that the actual absorption troughs 
comprise more than one intrinsic velocity component and are broader. In fact, Kaspi et al. 
(2001, 2002) measured directly v turb = 500 ± 170 km s" 1 . In fitting the spectrum, we assume 
a covering factor of unity, which is an approximation for the value of 0.8 obtained by Gabel 
et al. (2003), who fit for the covering factor of the X-ray observed components in the UV 
spectra. The present velocity-blended spectrum is not sensitive to this difference. Emission 
lines are modeled by Gaussians. None of these approximations hamper the good estimates 
we get for the total column density in each ion and for the total flux in each emission line, 
which are the physical quantities of interest here. 

The best-fit spectral model for the O region is given in Figs. 5 and 6, highlighting 
respectively the Ka and K/3 (and higher) transition regions. The good fit, including that for 
the high-order lines of the He-like series below 19 A, is apparent in the figures. Absorption 
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by all charge states from O VIII down to O IV is clearly detected. It is obvious from the 
figures that it is crucial to fit all charge states simultaneously, including both absorption and 
emission, in order to obtain the correct physical measurements. The various ions are found 
outflowing at (indistinguishable) centroid velocities of —470 — 800 km s _1 with respect to 
the systemic velocity of NGC 3783. The fitted parameters for each of these ions are given in 
Table 1. 

Evidently, the broad range of charge states observed along the line of sight is ubiquitous 
to Seyfert outflows (Sako et al. 2001; Behar & Netzer 2002; Blustin et al. 2003; Steenbrugge 
et al. 2003). It is also observed in Seyfert 2 galaxies (Kinkhabwala et al. 2002). If all 
charge states are present at similar distances from the central source, as observed directly in 
Seyfert 2's, it implies a strong density gradient of 3-4 orders of magnitude (e.g., for O IV to 
O VIII) over relatively short distances. In turn, comparable column densities together with 
a strong density gradient would imply that the highly-ionized gas occupies a much larger 
volume than the less-ionized gas. This picture is consistent with a wind that emanates from a 
dense object (e.g., accretion disk or bloated stars). It then accelerates, expands and becomes 
hotter and more ionized. Pressure equilibrium also favors the hotter gas being more tenuous 
and therefore more extended. 

In addition to the high charge states, weak absorption by O I and O II may also take 
place in (the rest frame of) NGC 3783 as indicated in Fig. 5. However, blending with other 
ions hampers an unambiguous detection. The O I Ka line in NGC 3783 overlaps with the 
strongest absorption lines of Ar XVI at 23.529 A (Lepson et al. 2003). The O II Ka line 
in NGC 3783 can not be resolved from the local galactic Ka line in our own galaxy. O III 
is not detected. Nevertheless, since we simultaneously fit many lines from each ion, we are 
able to overcome the problem of blending to some extent and to provide upper limits to the 
ionic column densities. Of course, these constrains are not as tight as those provided by UV 
measurements for the same ions, when those are available (Kraemer, Crenshaw, & Gabel 
2001). 

The errors on the outflow velocities are relatively large due to the wavelength channel 
widths of the RGS. Estimating the uncertainty to be bound by 3 channels implies an error 
range of ± 16.5 mA corresponding at 22 A to ± 230 km s _1 as quoted in Table 1. The 
best-fit column densities are also given in Table 1. The column densities of O VII and O VIII 
are constrained most tightly by their detected photoelectric edges. An additional constraint 
is imposed by fitting entire series of lines (see also Fig. 6). The values of (1.0 ± 0.3) and (4.0 
±1.0) x 10 18 cm -2 obtained for the column densities of O VII and O VIII, respectively, are 
in excellent agreement with the values of 1.1 and 4.3 x 10 18 cm -2 obtained by Kaspi et al. 
(2002). 
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The O VI ion is of particular interest for the possible connection between the UV and X- 
ray absorbers. Inner-shell O VI absorption can be observed in the present spectrum although 
only marginally. We detect the Ka line at a significance of about 2 sigma over a total of 
4 wavelength channels, but the errors on the data points are still too large to see the line 
profile clearly (compare data with continuum model in Fig. 5). The feature at 19.33 A can 
be associated with the O VI Kf3 line, although it is mostly due to O V K7 (see Fig. 6). 
The derived column density of O V is much higher than that of O VI. This is due to the 
closed shell structure of O V and to the pervasively low fractional abundance of the Li-like 
species O VI (< 0.25, e.g. Kallman & Bautista 2001). To that end, O VI and Li-like ions 
in general are perhaps the least favorable vestige for probing the ionized outflow in AGN. 
The precise restframe wavelengths of the Ka inner-shell transitions in O ions have yet to be 
measured in the laboratory and are still somewhat in debate (Behar & Kahn 2002). Here, 
we use wavelength values for O VI Ka and K(3 of 22.01 and 19.34 A, respectively. With a 
fitted turbulent velocity of 170 km s _1 , we obtain only an upper limit to the O VI column 
density of 2.0 x 10 16 cm~ 2 . This is loosely consistent with Gabel et al. (2003), who find a 
lower limit of 0.63 x 10 16 cm~ 2 for the total O VI column density in the two corresponding 
UV components. We note that there was a hint to the Ka line of O VI in the 900 ks HETGS 
spectrum (Kaspi et al. 2002), but there, the lower statistical quality did not allow for a 
definite identification. 

The present high-quality spectrum reveals an absorption feature at 21.64 A (21.43 A in 
Fig. 5 at the rest frame of NGC 3783). We identify it tentatively as absorption by O VII in 
the intergalactic medium (IGM). Given this interpretation, the gas would be at a redshift of 
z = 0.0018 ± 0.0008. The absorption parameters deduced for this line are given in Table 1. 
Another possible identification for this line could be Ca XVI (21.49 A in the rest frame) in 
the outflow of NGC 3783, but lines by neighboring charge states of Ca are not as obvious 
in the spectrum, which casts doubt on the identification of Ca. Since this feature may not 
be related to NGC 3783, we defer a more careful analysis of it to a separate paper. There 
is a hint of this line in the 900 ks HETGS spectrum, but there it seems to be at 21.60 A 
(as opposed to 21.43 A), where local O VII would be expected. In any event, due to the 
lower statistical quality of the HETGS observation in this spectral band, Kaspi et al. (2002) 
conclude that it is insignificant. 

2.3.2. Emission 

In addition to the numerous absorption lines, a few bright, narrow emission lines are 
present in the spectrum. In the O region of the spectrum, the brightest emission lines are 
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O VIII Lya at ~ 19 A and the Hea triplet of VII at ~ 22 A (the resonance, intercom- 
bination, and forbidden lines - r,i,Sz /). The radiative recombination continuum (RRC) of 
O VII is also detected, although its location at 16.8 A in the midst of the Fe-M 2p - 3d 
UTA makes a precise determination of its flux rather difficult. In fact, it is much better 
resolved with HETGS (Fig. 6 in Kaspi et al. 2002). The weaker O VII and O VIII K(3 as 
well as higher-order emission lines have been identified with HETGS (Kaspi et al. 2002), but 
here the lower spectral resolution impedes the measurement of their flux. The measurable O 
emission features are listed in Table 2. Fitting for a uniform velocity width for all emission 
lines, we obtain a turbulent velocity of 740 ± 140 km s" 1 (FWHM 1220 ± 230 km s" 1 ). Kaspi 
et al. (2002) measured a FWHM of 930 ± 350 km s^ 1 for the most isolated O VII /-line. 
Previous optical measurements by Appenzeller & Ostreicher (1988) of narrow, forbidden 
emission lines of highly- ionized species in NGC 3783 show a consistent linear increase in 
velocity width with ionization energy. The presently measured width of 1220 km s" 1 FWHM 
for O VII fits that trend perfectly (see their Fig. 13). For O VIII, we are not able to pro- 
vide an independent measurement of the emission line width, because there are no O VIII 
non-absorbed (forbidden) lines in the spectrum. 

As evident from Figs. 5 and 6, the emission/absorption model fits the data (including 
the intricate line profiles) remarkably well. This gives confidence in the individual (de- 
blended) line fluxes obtained. All of the measured emission lines are slightly redshifted by 
about 300 ± 230 km s -1 . Owing to the large errors, this is still consistent with the mean 
shift of 130 ± 290 km s^ 1 measured for the emission lines with HETGS (Kaspi et al. 2002). 
These redshifts together with the complex P-Cygni type profiles seen for the resonance lines 
are suggestive of a bipolar outflow as discussed in detail below. Also presented in Table 2 
are the line fluxes measured with HETGS. The current values are in good agreement with 
the HETGS values. The largest discrepancies are found for the resonance lines, which have 
complicated P-Cygni type profiles [Lya and Hea(r)]. This is due to the different methods 
used for measuring the line fluxes. For example, Kaspi et al. (2002) did not fit for overlapping 
absorption and did not require the widths of the three O VII lines to be uniform. Nonetheless, 
the general agreement between the two measurements implies no time variation as expected 
for the narrow line region (NLR). 

3. Discussion 

3.1. Relation Between the X-Ray Emitter and Absorber 

To summarize the X-ray outflow kinematics in the rest frame of the AGN as mani- 
fested by the absorption and emission lines: The emission lines are redshifted by about 
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300 ± 230 km s" 1 and have FWHM of 1220 ± 230 km s" 1 (present work). The absorption 
lines are blueshifted by about -600 ± 150 km s" 1 and have a FWHM of 820 ± 280 km s" 1 
(Kaspi et al. 2002). Although it is impossible to constrain the exact shape of the X-ray line 
profiles and in particular their broad wings, due to the crowded nature of this waveband 
and the overlapping emission and absorption, it is obvious that high-velocity gas (up to 
~1000 km s" 1 ) is both emitting and absorbing X-rays. For the resonance lines, P-Cygni 
type profiles are observed (Figs. 5 & 6). These line profiles are characteristic of a bipo- 
lar, non-collimated outflow in conjunction with absorption by the outflow along the line of 
sight. It is important to point out that we have only approximated these profiles by emission 
and (slightly offset) absorption lines, but the comparable velocity widths obtained for the 
absorption and emission lines already supports the biconical interpretation. 

We wish to further compare the X-ray line emission and absorption in NGC 3783. In 
the context of Seyfert 2 galaxies, it has been shown that intensity ratios among the He-like 
emission lines can be used to decisively determine the column density through the emitting 
plasma (Behar et al. 2001; Kinkhabwala et al. 2002) even when the absorption spectrum is 
not available to us directly due to orientation & obscuration effects. The present emission 
line fluxes give ratios of f/r — 1.5 ± 0.3 and i/r = 0.5 ± 0.2. With the velocity width of 
740 km s _1 , this f/r value is obtained in our model at an O VII ionic column density of (1.1 
± 0.35) x 10 18 cm" 2 . The i/r ratio provides just a lower limit of 1.3 x 10 18 cm" 2 . These 
ranges are both in good agreement with the ionic column density of (1.0 ± 0.3) x 10 18 cm -2 
measured directly from the absorption lines and edge of O VII. 

Next, it would have been useful to compare the fractional ionic abundances derived from 
the emission lines with those derived from the absorption. However, direct comparison is 
not available since recombination-driven lines are due to the next ionized state (e.g., O IX 
for Lya). Hence, the O VII & O VIII emission line intensities reflect the novni/noix 
abundance ratio, while the absorption column densities reflect riovn l n ovin- Nevertheless, 
the comparable emission line fluxes and column densities of O VII and O VIII suggest that 
all highly-ionized states are fairly abundant in the emitting plasma, just as they are in the 
absorbing plasma. Moreover, it is possible to find an ionization parameter £ = 30 erg s" 1 cm 
that reproduces just the right amount of O VII and O VIII emission and absorption, which 
we observe. Clearly, we are seeing a very broad distribution in £, but for simplicity, this 
single value can be used to describe the O VII & O VIII spectral features. In summary, the 
velocity widths, the column densities, and the fractional abundances all make it tempting to 
interpret the X-ray emitter and absorber as arising from the same gas. 
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3.2. Location and Geometry of the X-Ray Plasma 



The measured emission line flux Fji (e.g., of an ion) can be related to the volume (V) 
emission measure EM as follows: 



The computed line power Pji is defined through the line emissivity n e n q Pji (i.e., number 
of line photons emitted in the plasma per unit time and per unit volume). Here, n q is the 
number density of ions with charge q. The line power takes into account all line-driving 
mechanisms and in particular recombination, radiative cascades, and photoexcitation. The 
fraction of ions with charge q (n q /Hn q ') is denoted by f q . For recombination-driven emission 
lines such as the present ones, the most relevant fraction is f q +i- In eq. (1), A Q represents 
the oxygen elemental abundance relative to Hydrogen. The distance to the source d is taken, 
for the case of NGC 3783, to be 41.8 Mpc. In a fully ionized cosmic plasma, we can also 
relate the electron density to the Hydrogen density by n e = 1.2 uh- On the other hand, the 
column density N H can be written as follows: 



The measured ionic column densities iVj are related to Nh by Ni= f q AoNu. The minimal 
and maximal distances for absorbing gas in the outflow are denoted by r min and r max . To 
obtain the right hand side of eq. (2), we used n e = L/(£r 2 ). 

We focus here on the O VII and O VIII spectral features as the tracers of the very 
highly-ionized gas simply because they exhibit the best measurable features in our spectrum. 
However, since Seyfert 2 observations show that all charge states co-exist and overlap in space 
along the outflow, we expect the present analysis of O ions to be applicable to the entire 
highly-ionized plasma. We choose an ionization parameter of £ = 30 erg s" 1 cm, where 
ionization balance calculations give novm l n oix — 1-7 (in accord with the emission line 
ratio) and novm / n ovn — 4 (in accord with the absorption column densities). This implies 
f q values of 0.14, 0.57, and 0.29, respectively for O VII, O VIII, & O IX. Again, a single 
value of £ is assumed in this section just for simplicity. The observed broad distribution in 
£ makes it clear that each ion exists mostly around its typical range of £. 

Equations (1) & (2) can not be solved for the most general case and some approximations 
need to be made. In the following we explore the two extreme scenarios: r max = r min and 




(1) 





r. 



max 



mm- 



3.2.1. Localized Scenario 



First, we explore the possibility of a well localized X-ray plasma, i.e., one in which the 
typical size Ar = r max — r m i n is much smaller than the average distance r from the ionizing 
source. In that case, the localization in r immediately imposes a narrow distribution in n e 
(through £). Eqs. (1) and (2) can then be rewritten simply as: EM = nnn e Ar 3 and N H = 
(n H /n e )(L I £)(Ar /r 2 ) . Even now, these equations can not yet be solved separately. However, 
if one identifies the emitter with the absorber, both equations can be solved simultaneously. 
Indeed, by using eqs. (1), (2), and the definition of £ to solve for r, Ar, and n e one would 
expect to obtain a rough idea of these quantities. However, using the present measurements 
for the O VII forbidden line and for the O VIII Lya line, in both cases, we obtain values 
of r = Ar = 10 pc. This suggests that the confined geometry (Ar <C r) is inadequate 
for describing the highly-ionized outflow of NGC 3783, at least as long as the line emission 
and absorption are due to the same outflow constituent. All of the measured and derived 
quantities are summarized in Table 3. 



3.2.2. Ionization Cone Scenario 

Inspired by the Chandra X-ray images of Seyfert 2 galaxies (Sako et al. 2000; Young, 
Wilson, & Shopbell 2001; Brinkman et al. 2002; Ogle et al. 2003), we propose an alternative 
geometry of two ionization cones with an opening solid angle Q. In NGC 1068, AGN driven 
X-rays are detected as far as 500 kpc away from the nucleus (Brinkman et al. 2002). Further- 
more, Young et al. (2001) & Ogle et al. (2003) show that the X-ray NLR in NGC 1068 overlaps 
on the plane of the sky with the optical NLR. Extended ionization cones in Seyfert l's, only 
oriented along the line of sight, would substantiate the unified theory of AGN (Antonucci 
& Miller 1985) in the X-ray regime. Unfortunately, it is impossible to distinguish between 
clumpy clouds and ionization cones in NGC 3783 by means of X-ray imaging, as even the 
best angular resolution available with Chandra (0".5) corresponds to 100 pc at the distance 
of NGC 3783. Therefore, we have to make use of the spectrum. 

In the ionization cones scenario, the plasma covers a large range of distances (and 
densities) and it is necessary to revert to the integral forms of EM and in eqs. (1) and 
(2) with dV = 2Vtr 2 dr. Note that even within the framework of the extended geometry and 
even when r max ^> r min , for a given £ most of the absorption occurs close to r min , while the 
actual size of the cones (i.e., r max ) can not be well constrained by the observation. Indeed, 
eqs. (1) and (2) can now be rewritten as: 
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n H L 



(3) 



min — 



mm 



max 



and 



1 EM £ 

2 A^"Z 



(4) 



In this scenario, since eq. (3) is attained from the column density alone without the 
emission measure, the result for r min is independent of the emission. Thus, it is independent 
of the assumption that the X-ray emitter and absorber represent the same plasma. Eqs. (3) 
& (4) demonstrate the dependence of the derived values on the assumption for £. Thus, 
if O VII actually forms at £- values somewhat lower than we assume (30 erg s _1 cm), we 
would be slightly underestimating r min and overestimating Q, and vice versa. We stress 
that the focus on a certain ion or on a specific value of £ constitutes a complementary and 
totally consistent approach with the observation of a broad range of charge states. Again, 
analogous to the Seyfert 2 case, each charge state could be present throughout the cone, only 
at different densities, while at the same time different charge states are co-located thanks to 
the steep density gradients at each point in the cone. The X-ray observation simply traces 
O VII and O VIII wherever they are along the cone. 

We now substitute the measured values into eqs. (3) & (4) using the data for the O VII 
forbidden line and for the O VIII Lya line. These values as well as the results for r min and 
for Q are summarized in Table 3. The electron density is again n e = L/(£r 2 ), but here it 
decreases rapidly (oc r~ 2 ) along the cone. In Table 3 we give the maximal n e value attained 
at r min . It can be seen that the results for both O VII and for O VIII suggest consistently 
that the X-ray emitter/absorber lies at r min = 10 pc. The maximal densities are of the 
order of 600 cm -3 . These results are not very different from those obtained by Shields & 
Hamann (1997), who do not assume a cone geometry, but give values of n e > 50 cm" 3 and 
r m in < 30 pc. 

The typical opening angle Q/2it we get for each of the two cones is about 15 - 25 %, 
fully consistent with the fraction of Seyfert galaxies that show a highly-ionized absorbing 
outflow (~50% of the Seyfert l's). The presently proposed geometry with the high r min and 
Q values makes for a large scale outflow that from our perspective would entirely blanket the 
central AGN region including the BLR, provided that the ionization cone is uniformly filled. 
Since Gabel et al. (2003) find that at some radial velocities the BLR is only partially covered 
by the UV outflow, identifying the UV absorber with the extended X-ray ionization cone 
would require that some of the BLR be observed directly, i.e., through (velocity-specific) 
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holes in the cone. We note that Steenbrugge et al. (2003) estimate for NGC 5548 much 
smaller opening angles, Q/2tt < 0.5%. However, that result was obtained by assuming a 
mass loss rate that is equal to the mass accretion rate (with a luminosity efficiency of 5%). 
If this assumption is relaxed, extended ionization cones may be applicable for NGC 5548 as 
well. Indeed, the high values for Q in NGC 3783 imply that the outflow mass loss rate exceed 
the accretion rate as discussed in the following section. 



3.2.3. Mass Estimates 

By integrating the mass density fxnuvnu (where fx = 1.3 is the average atomic weight 
for a solar- abundance plasma) over the ionization cone (or clouds), it is possible to estimate 
the total mass M and similarly the mass loss rate M in the outflow. These quantities 
are listed in Table 3 for both the localized and ionization cone models. In the ionization 
cone model, the total mass depends on r max — r min , which is not well constrained by the 
observation. Therefore, the mass in Table 3 is normalized to r max = r max — r min . Taking 
an outflow velocity of —600 km s -1 , we obtain a total mass loss rate of up to a few solar 
masses per year. We compare this with an estimate of the mass accretion rate. With a 
bolometric luminosity of 4.5 xlO 44 erg s _1 (Markowitz et al. 2003) and a gravitation-to-light 
conversion efficiency of Lt> i/(M acc c 2 ) = 5%, NGC 3783 would be accreting about 0.2 solar 
masses per year. If all these numbers are correct, then NGC 3783 is accreting only about 
10% of the mass it is ejecting. Although this is a surprisingly low accretion efficiency, we are 
not aware of any hard evidence against it. Furthermore, ultra-high velocity outflows (~0.1c) 
have been discovered recently in several AGN. It has been postulated that the high mass 
loss rate observed in these sources is due to their high accretion rate L ho i/ L Eddington (Pounds 
et al. 2003). Narrow line Seyfert l's such as NGC 3783 could be the low-mass facet of this 
phenomenon. 

It is worth stressing that the errors given in Table 3 reflect only those due to the presently 
measured iVj & values. Indeed, the present mass estimates should be treated with more 
caution as the X-ray outflow of NGC 3783 is undoubtedly more complex than a simple, 
continuous cone, which was assumed here. Nonetheless, these estimates are good enough 
to rule out radically different scenarios, such as the highly-localized one. Most likely, the 
true picture for NGC 3783 is somewhat less extreme than r max ^> r min and therefore r min 
could be somewhat less than 10 pc, but certainly not much less, given the timing constrains 
obtained in Sec. 2.2 of r min > 0.5 pc for O-K and r min > 2.8 pc for Fe-M. 
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3.2.4- Comparison With Other Observations 

The extended ionization cones stand in contrast with interpretations of two other 
Seyfert 1 galaxies. Netzer et al. (2002), based on observed recombination times, place the 
ionized absorber in NGC 3516 about 0.2 pc away from the nucleus. Kriss et al. (2003) and 
Blustin et al. (2003) even suggest that part of the ionized absorber in NGC 7469 is closer 
to the nucleus than the BLR. In both of these cases, the observers find that the UV and 
X-ray absorbers are connected. More importantly, the present results are in stark contrast 
with the parameters derived by Gabel et al. (2003) for NGC 3783 based on UV absorption, 
which according to those authors occurs at sub-parsec distances in tiny clouds the size of 
Ar = 10 10 cm. At the heart of the argument of Gabel et al. (2003) lies the absorption by 
the C III multiplet at 1175 A detected in the high- velocity UV component (-1300 km s _1 ). 
Following Bromage et al. (1985), Gabel et al. (2003) take this feature as evidence for high- 
density gas (> 10 9 cm -3 ) in this component, which inevitably restricts it to clumps very 
close to the source. We wish to point out that the 2s2p triplet complex from which those 
C III lines are absorbed, can in fact be populated at lower densities. Only one fine-structure 
level in the triplet requires high densities (see, e.g., Fig. 1 in Bhatia & Kastner 1993), but 
absorption from the other two levels of the triplet may be possible at low densities and might 
be producing the C III lines at 1175 A observed by Gabel et al. (2003). 

On the other hand, it may also be possible that dense clouds reside at the base of the 
ionization cone where the outflow is clumpy, less ionized, and thus absorbs UV light. Further 
out in the outflow, perhaps, the ejecta evaporate, become much more tenuous and conse- 
quently more highly-ionized. These more extended parts of the outflow might then constitute 
the X-ray absorbing/emitting gas. However, following Kaspi et al. (2002) and Gabel et al. 
(2003), we note that the very similar velocity structure in the X-ray and UV absorption sys- 
tems strongly suggests that the two coincide and most likely are two different manifestations 
of the same outflow. Finally, based on UV variability and recombination times, Shields & 
Hamann (1997) also arrive at the conclusion that the X-ray and UV absorbers are more than 
10 pc away from the ionizing source. 

A paramount difference between the UV-derived and X-ray-derived parameters of the 
outflow lies in the corresponding M and M values. In the Gabel et al. (2003) picture, the 
calculated M and M values are more than 10 orders of magnitude smaller than the values 
derived here from the X-ray measurement (see Table 3). According to the UV interpretation, 
the outflow carries insignificant mass and energy out of the AGN system, while according 
to the present X-ray analysis, the outflow must play a major role in the mass and energy 
budgets of the AGN. 
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4. Conclusions 

The present 280 ks observation of NGC 3783 leads to a few interesting conclusions, which 
can be summarized as follows: 

1. Despite the continuum increasing by up to a factor of 2 during the 280 ks observation, 
no significant effect on the absorption structure has been detected. The fact that the 
absorbing plasma did not respond to the increase in flux within, say, 1 day places it at 
least 0.5 pc (O-K) to 2.8 pc (Fe-M) away from the continuum source, regardless of its 
geometry or physical situation (e.g., density). 

2. The interpretation of the X-ray line absorbing and emitting plasma being closely con- 
nected, or actually being manifestations of the same outflow, is very tempting given the 
P-Cygni type profiles as well as the similar velocity widths, inferred column densities, 
and ionization distributions. 

3. Assuming that the two components are connected rules out the localized geometry 
(Ar <C r and n e = 10 9 cm -3 ) for the X-ray plasma. This geometry was suggested by 
Gabel et al. (2003) for the UV absorbing outflow. 

4. Instead, we propose diffuse gas in the form of ionization cones. These cones may not 
be closer than about 10 pc from the ionizing source and must therefore have much 
lower densities (a few 100 cur 3 and decreasing with distance). In this case, we find 
that the mass loss rate in the outflow is significant, of the order of a solar mass per 
year. 

5. The conflict with the clumpy cloudlet description of Gabel et al. (2003) could possibly 
be due to erroneous use of the C III multiplet at 1175 A as a density diagnostic, but 
it could also be resolved if the UV absorber lies at the base of the ionization cone in 
the form of clouds. This constituent then represents a transient phase of the gas as 
it expands, ionizes, and transforms into the much more extended and more tenuous 
X-ray plasma. 

This work is based on observations obtained with XMM-Newton an ESA science mis- 
sion with instruments and contributions directly funded by ESA Member States and the 
USA (NASA). EB acknowledges helpful discussions on AGN phenomenology with Ari Laor. 
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Table 1. 


Ions seen in absorption in 


the RGS spectrum 


of NGC 3783. 


Ion 


Outflow Velocity [km s _1 ] 


vturb [km s^ 1 ] a 


Ni [10 18 cm" 2 ] 


VIII 


-550 ± 230 


170 ± 50 


4.0 ± 1.0 


VII 


-660 ± 230 


170 ± 50 


1.0 ± 0.3 


VI 


-800 ± 230 


170 ± 50 


0.01 ± 0.01 


V 


-470 ± 230 


170 ± 50 


0.18 ± 0.03 


IV 


-770 ± 230 


170 ± 50 


0.08 ± 0.02 


III 






< 0.005 


II b 


± 230 


A;T;=0.1 eV c 


0.004 ± 0.01 


I d 


± 230 


A;T i= 0.02 eV 


0.025 ± 0.035 


VII Gal. c 


z = 0.0018 ± 0.0008 


40; kTi=0.1 keV 


0.084 ± 0.03 



a Turbulent velocities of O IV - O VIII were tied together in the fit. 
b Blended with O I in our galaxy. 
c Ion temperature. 

d Blended with Ar XVI in NGC 3783. 
c Tentative identification. 
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Table 2. Oxygen emission lines in the RGS spectrum of NGC 3783. 



Ion 


Observed 


Velocity Shift 


Velocity Width 


Line Flux 


Line Flux 




Wavelength 


[km s -1 ] a 


[km s" 1 ] b 


[10~ 5 cts s" 1 cm" 2 ] 


[10~ 5 cts s' 1 cm" 2 ] 




[A] 








(Kaspi et al. 2002) 


VIII Lya 


19.173 


330 ± 230 


740 ± 140 


4.5 ± 0.8 


2.32 ± 0.74 


VII r 


21.833 


290 ± 230 


740 ± 140 


7.8 ± 1.5 


4.58 ± 2.38 


VII % 


22.034 


290 ± 230 


740 ± 140 


4.1 ± 1.3 


3.52 ± 1.82 


VII / 


22.335 


290 ± 230 


740 ± 140 


11.5 ± 1.6 


10.13 ± 2.74 



a In the system of NGC 3783 and tied for all O VII lines. 
b Tied for all lines. 
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Table 3. Measured and derived parameters for the X-ray emitter /absorber in NGC 3783 

as traced by VII & VIII. 



Parameter 



O VII (/) O VIII (Lya) UV (Gabel et al. 2003) 



s 1 cml 



f [erg 

L [10 43 erg s 
Ni [10 18 cm" 2 

u 

fq+1 

Ao 

N H [10 22 cm 
d [Mpc] 

Fji [10~ 5 cts s^ 1 cm 
Pji [10~ 12 s" 1 cm 3 
EM [10 64 cm" 3 ] 



-21 



-21 



30 
1.5 
1.0 ± 0.3 
0.14 
0.57 
5xl0~ 4 

I. 43 ± 0.43 

41.8 

II. 5 ± 1.6 

3.9 
2.2 ± 0.3 



30 
1.5 
4.0 ± 1.0 
0.57 
0.29 
5xl0" 4 
1.40 ± 0.35 

41.8 
4.5 ± 0.8 

5.7 
1.2 ± 0.3 



r [pc] 
Ar [pc] 
n e [cm- 



Localized Clouds Scenario (r max — r mi , 
16.3 ± 7.3 12.3 ± 4.9 
9.9 ±3.0 8.1 ±2.1 
200 ± 180 340 ± 270 



= Ar < r) 

<0.26 
3.3 xlO- 9 
10 9 



v out [km s" 1 ] 
M [10 37 g] 
M [10 33 g yr- 1 ] 
n = (Ar/r) 2 [str] 



-600 ± 150 -600 ± 150 
1.0 ± 1.3 1.0 ± 1.1 
0.6 ±0.7 0.8 ±0.8 



-550 - -1300 
1.8 xlO" 22 
3.4 xl0~ 13 
1.6 xlO" 16 



[PC] 



n e (r min ) [cm 3 ] 



Q/271 



Ionization Cones Scenario (r max r ri 
9.4 ± 2.8 9.6 ± 2.4 
590 ± 350 560 ± 280 
0.24 ± 0.08 0.13 ± 0.05 



<0.26 
10 9 
2.6 xlO" 17 



v out [km s 4 ] 
M/r max [10 36 g pc- 1 ] 
M [10 33 g yr- 1 ] 



-600 ± 150 
8.4 ± 2.6 
5.2 ± 2.0 



-600 ± 150 
4.7 ± 1.8 
2.7 ± 1.0 



-550 - -1300 
M = 1.8 xlO 15 g 
3.4xl0- 13 
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Fig. 1. — The RGS light curve (RGS 1 and 2 combined) during the present 280 ks observation 
in time bins of approximately 10 minutes. Different segments of the observations have been 
labeled. 
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Fig. 2. — The RGS spectra during the four time segments defined in Fig. 1. The data are 
plotted in bins of about 0.1 A to accentuate broad structure. Hence also the small error 
bars. The lines between the data points are plotted just to guide the eye. Apart from the 
different flux levels, the shape of the absorption structure appears to be remarkably similar 
in all four spectra. In particular, the broad absorption features due to the Fe-M 2p-3d UTA 
(16 -17 A) and due to the O-K edges (14.2 & 16.8 A) do not appear to vary. The quality of 
the data is not good enough to enable detection of variability in the narrow lines. 
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Fig. 3. — Ratio of the "high" phase and "low" phase spectra rebinned to about 0.1 A. 
No appreciable broad deviations from the mean (horizontal line) are found implying no 
significant change in the photoelectric edges (marked above the data). The gaps in the 
plot are instrumental. Note that even with no change in absorption, the ratio around strong 
emission and absorption lines (e.g., at 13.5, 14.2, 15, and 19 A) is not expected to be perfectly 
flat, since deep absorption troughs suffer from low S/N and emission lines do not cancel out 
in the ratio. 



-26- 




Fig. 4.— The total 280 ks RGS spectrum of NGC 3783. The slab model from Blustin 
et al. (2002), which was fitted to the first 40 ks RGS observation, is over-plotted with 
no adjustments, including no re-normalization. The over all agreement is very good, but 
noticeable discrepancies are seen in the details of the O-K region between 14 and 24 A. 
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Fig. 5. — Spectrum of NGC 3783 in the region of the Ka lines of O I - O VII as obtained 
with a 280 ks RGS observation. The thick curve represents the best-fit model and the thin 
curve is the continuum model only with galactic and instrumental absorption. 
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Fig. 6. — Spectrum of NGC 3783 showing the region of O VIII Ka as well as higher order 
lines (K/3, K/y, etc.) of lower charge states of O as obtained with a 280 ks RGS observation. 
The thick curve represents the best-fit model and the thin curve is the continuum model 
only with galactic and instrumental absorption. 



